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Deviations from Fermi liquid behavior are well documented in the normal state of the cuprate 
superconductors, and some of these differences are possibly related to pre-formed pairs appearing 
at temperatures above T c . In order to test these ideas we have investigated the attractive Hubbard 
model within a self-consistent, conserving ladder approximation. In this version of the theory, no 
feature is present which can be related to the pseudo gap found in the high-T c materials. Further, the 
interactions between two-particle bound states change the physics of the superconducting instability 
in a profound fashion, and lead to a completely different phenomenology that one predicts based on 
the non-self-consistent version of the same theory. 
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Numerous recent experiments (e.g., ARPES, optical, 
neutron scattering, NMR) have shown that in the high 
T c cuprate superconductors a so-called pseudo-gap is 
present jj] ; this is further supported by the extensive the- 
oretical work of Randeria and collaborators 0] . This has 
led to proposals that electron pairs form at temperatures 
well above the superconducting transition temperature. 
However, possibly due to phase fluctuations, a macro- 
scopic phase coherent wave function is not formed, so 
superconductivity is not encountered until lower temper- 
atures ||. 

This physics motivates our study of the attractive Hub- 
bard model. The Hamiltonian for this system is 

H = ~tJ2 ( c U c J>+ h - c -) - \ U \ I>*,t n U (!) 

where the lattice sites of a 2-D square lattice are labeled 
by {i}, the lattice fermion operators are denoted by c% i<T , 
and neighbouring sites are represented by (if). For any 
nonzero | U |, two-particle bound states appear, and are 
physically related to a pair of electrons lowering the sys- 
tem's energy when they exist on the same lattice site. 
This is certainly the simplest example of a model Hamil- 
tonian which allows for one to study the interactions be- 
tween such electron pairs. 

We employ the Brueckner-Hartree-Fock theory (or lad- 
der approximation) which neglects all crossing diagrams 
and which is therefore valid only at low electronic den- 
sities. A preformed pair manifests itself in this theory 
as a two-particle bound state. In the non-self-consistent 
version of the theory such a pair of electrons always has 
an infinite lifetime. 

By analyzing a fully self-consistent formulation of this 
theory, we are able to investigate preformed pairs with 
arbitrary lifetime. That is, if one solves this problem in 
this approximation, one includes pair-pair interactions. 
Haussmann has argued |i| (when he investigated a three- 
dimensional continuum model with the same physics) 
that the interactions between pairs will be, in first or- 
der, of a repulsive nature. This leads to the dissociation 
of the pairs causing its finite lifetime. 

The equations for the Green's function (G), pair sus- 
ceptibility (x), four-leg vertex function (r), and self en- 
ergy (E), in a conserving approximation are well known: 

G(k,uo n ) = (G (k,^„)- 1 -E(k,^„)) _1 (2) 
X(K, iO n ) = — ^2 G ( K ~ k ' ^ - iu m )G(k, iuo m ) (3) 

m.k 

r(K, iQ n ) = | U | / (1+ | U | X (K, «„)) (4) 
S(k, iu> n ) = — ^2 r ( k + Q= iuj n + iw„)G(q, iu m ) (5) 

m,q 

where the wave vectors and Matsubara frequencies have 
their usual meaning. In a non self-consistent theory, 



one replaces the full Green's functions G in Eqs. (3,5) 
with the noninteracting Green's functions Go- The usual 
Thouless criterion begins with the non self-consistent the- 
ory and identifies an instability Q of the normal state 
in terms of poles of the pair susceptibility (in a weakly 
interacting 3-D system, the Thouless criterion becomes 
identical to the equation which determines T c in the BCS 
theory). The changes in this phenomenology that arise 
when one includes pair-pair interactions may be an im- 
portant component of the high T c problem. 

For full self consistency, this set of equations has to 
be solved iteratively. Since such solutions are difficult 
and time consuming to obtain, we have investigated a 
simple approximation that allows for extensive numeri- 
cal investigations of the resulting equations. To be spe- 
cific, during the first step of the iteration process leading 
to self consistency, we make an approximation for the 
pair susceptibility as being equal to the k-average (de- 
noted from now on as over-lined quantities, e.g. T) of 
the noninteracting pair susceptibility. We only calculate 
the k-averaged pair susceptibility during subsequent it- 
erations to self consistency. This leads to the following 
expressions 

F(ifin) = 1 ]T r(K, in„) «| u | /(i+ 1 u | (6) 

K 

G(k,iu n ) « (G (k, i ^„)- 1 - E(io; n )) _1 (7) 

Besides the computational simplicity of this approach, 
one important additional advantage of this technique 
arises from the ease with which low temperatures may be 
reached for two and higher dimensional systems. Other 
methods have great difficulty reaching the low tempera- 
tures for which superconducting properties may become 
apparent, potentially making our approach quite valu- 
able. 

Our k-averaged approximation becomes accurate in 
any of the following limits: Quite trivially, (I) At large 
temperatures all correlations are thermally washed out, 
and for small bandwidths one obtains the atomic limit. 
In both cases, k-dispersion of the pair susceptibility is 
irrelevant. (II) For large U, T is basically determined by 
the pole of the bound state at \ — 1/f which means 
that a weakly dispersive two particle bound state is well 
separated from the continuum. So, in this case the av- 
erage over k-space is a good approximation, since the 
bound state below the continuum is now replaced by it's 
k-average. (Ill) Probably the most interesting limit in 
which our approximation is valid is that of large spa- 
tial dimensions. It was argued by Vollhardt and collab- 
orators (6) that for a system with large dimensions the 
k-dispcrsion of the self-energy vanishes. Such an approx- 
imation gives several results which are also valid in lower 
dimensions. This limit gives us a tool to compare our 
final result with well known results obtained from the 
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Kondo-impurity problem. This is due to the possibility 
of mapping the attractive Hubbard model, via a particle- 
hole transformation, onto a positive U Hubbard model in 
a magnetic field. The repulsive Hubbard model, in the 
limit of large spatial dimensions, becomes equivalent to 
the Kondo impurity problem 0. We will discuss more 
fully the regimes of validity of this approximation in a 
future publication. 

We begin by showing that for a one dimensional system 
the k-average method correctly reproduces the physics 
found in a fully self-consistent, conserving theory even in 
the intermediate coupling regime. In Fig. |l| we have plot- 
ted the variation of electron density with temperature 
for | U | being equal to the bandwith for \i = — t (note 
that the free electron band extends from — 2t to 2t). This 
quantity was calculated firstly for the self-consistent, con- 
serving theory obtained from the k-average method, and 
then from the full theory where k-dispersion was taken 
into account. The latter work involved 150 points along 
the imaginary frequency axis and 40 real-space lattice 
points — we found that these numbers gave fully con- 
verged results. In both cases the electron density de- 
creases monotonically with temperature, and there is no 
sign of a divergence (or strong increase) of the particle 
number. 

We note that this behaviour is in direct contrast to 
results obtained from non self-consistent calculations. In 
one and two dimensions, for low temperatures, the sys- 
tem is always unstable towards a condensation into an 
infinite lifetime two-particle bound state. For two di- 
mensions this was demonstrated for a non self-consistent, 
non conserving theory by Schmitt-Rink, et al. ||, and 
we Q have recently shown similar physics is encountered 
when a fully conserving theory is used. However, the 
self-consistent theories, at least to the minimum temper- 
atures that we have been able to access, do not show 
similar physics. 

In this brief report we want to focus on only one as- 
pect of the numerical results which we have generated 
from our k-averaged treatment of the self-consistent, con- 
serving treatment of the attractive Hubbard model. In 
Fig. H we show the density of states A(to) at different 
temperatures for a two-dimensional square lattice, again 
taking U | to be equal to the bandwidth, at a low elec- 
tron density n = 0.3 (n = 0.5 corresponds to half filling). 
We have also shown the results that follow from other 
approximations. 

The non self-consistent conserving approximation 
leads to a large fraction of the density of states being 
shifted to a peak below the one particle continuum - 
this peak is due to the infinite lifetime two-particle bound 
state. Between this peak and the one particle continuum 
lies a large gap, and this gap has been suggested by some 
to be the origin of the pseudo gap. 

However, our self-consistent, conserving, k-averaged 



calculation shows that the gap mentioned above disap- 
pears entirely, and that the only remnant from the two 
particle bound state is a flat, strongly lifetime broadened 
shoulder below the continuum. Further, instead of a gap, 
the density of states becomes a small maximum at the 
chemical potential, and with decreasing temperature this 
maximum increases. The existence of such a maximum 
can be understood via the analogy of the attractive Hub- 
bard model with the repulsive Hubbard model which in 
large spatial dimensions can be mapped onto a Kondo 
impurity model as mentioned above. In such problems a 
local maximum of the one particle density of states also 
occurs, as seen, e.g., by Georges et al. JTo| , who investi- 
gated the repulsive Hubbard model by using an iterated 
perturbation theory. 

We note that recent work for a d-wave pairing inter- 
action []llf has found similar results, but has been inter- 
preted quite differently. 

In conclusion, we have demonstrated the possibility of 
solving the attractive Hubbard model in two dimensions 
in the self-consistent, conserving ladder approximation 
by performing a numerically advantageous k-average. 
Due to self-consistency, we include pair-pair interactions 
which hinder the formation of an infinite lifetime bound 
state. In this way, the Thouless criterion is altered. Fur- 
ther, evidence for a pseudo-gap in the density of states is 
not obvious. 

We thank Frank Marsiglio for helpful comments. This 
work was supported by the DFG (Deutsche Forschungs- 
gemeinschaft), and the NSERC of Canada. 
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FIG. 1. A demonstration of the success of the k-averaged 
method for a one-dimensional system. Here, | U \ /t = 4, 
namely the Hubbard energy is equal to the bandwidth, and 
the chemical potential is fixed to be n/t = — 1. Then, the elec- 
tron density n is calculated as a function of temperature. The 
solid line shows the "exact" calculation for the self-consistent, 
conserving approximation, and the k-averaged result is shown 
by the dashed line. All energies are given in units of t. 
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FIG. 2. The density of states A(U>) that results for the 
two-dimensional system with n=0.3 and | U \ /t = 8 for 
four different temperatures. Note that the density of states 
develops a maximum at the chemical potential when the tem- 
perature is decreased. 
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